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We propose a family of master equations for local quantum dissipation. The 
master equations are constructed in the form of Lindblad generators, with the con- 
straints that the dissipation be strictly linear (i.e. ohmic), isotropic and translation- 
ally invariant. The resulting master equations are given in both the Schrodinger and 
Heisenberg forms. We obtain fluctuation-dissipation relations, and discuss the relax- 
ation of average kinetic energy to effective thermal equilibrium values. We compare 
our results for one dimension to the Dekker master equation [H. Dekker, Phys. Rep. 
80, 1 (1981)], which can be interpreted as a low length scale approximation of our 
model, as well as the Caldeira- Leggett master equation [A. O. Caldeira and A. J. 
Leggett, Physica (Utrecht) 121 A, 587 (1983)]. 
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I. INTRODUCTION 



Dissipative quantum systems and quantum brownian motion have long been of interest 
and a wide variety of techniques have been used to investigate them. We will not 

attempt a comprehensive review, but extensive reviews can be found in Refs. |IJ and 0. For 
linear dissipation, one wishes to obtain dynamics for a quantum mechanical system which 
are in a sense generalizations of the classical equations of motion for linear dissipation: 

mx = p, 

P = --P + F. (1.1) 

m 

The classical Langevin equations (and nonlinear generalizations) can be obtained by con- 
sidering a particle interacting with an oscillator bath. || Quantum mechanical calculations 
often begin with the same type of environment. For example, Caldeira and Leggett [|J use 
the influence functional technique of Feynman and Vernon || for a particle linearly coupled 
to a thermal bath of harmonic oscillators. The resulting (high temperature) master equation 
is given by 

An alternative approach by Dekker [[TJ, uses the quantization of complex classical canonical 
coordinates and quantal noise operators, resulting in the master equation: 

{p, p}] + ^(D xp + D px )[x, [p, p}} 

-^\p>\p>p]]--j!jr[ x >[ x >p]]- ( L3 ) 

Neither the Caldeira-Leggett nor the Dekker master equation is translationally invariant, 
and so cannot represent the effect on a quantum particle of a local interaction with a spa- 
tially homogenous environment. Caldeira and Leggett have pointed out that it is only 
strictly linear coupling (i.e. bilinear in system and oscillator bath coordinates) that can 



produce strictly linear dissipation. However, it may be true that other forms of coupling 
produce linear dissipation as an approximation which is valid over a wide range of conditions. 
Moreover, the dynamics of quantum coherence could be very different in a nonlinear (local 
coupling) versus a linear model ||[7| . It is the intent of this paper to show one method of 
constructing a master equation which reflects (in translational invariance) the symmetries 
of the environment and environment- system coupling which presumably give rise to the 
dissipative dynamics. 

The approach to constructing the master equation is essentially algebraic. We begin with 
a master equation of the Lindblad form: 

d P T\ i 

m = L[p] 

= ^-[H,p]+AL[ P ], (1.4) 

which, for bounded H, {V^}, generates a completely positive dynamical semigroup. || The 
generator for the corresponding Heisenberg picture is given by: 

dt dt 1 J 

= ~[H,0] + AL m [0]. (1.5) 

Following an approach used by Sandulescu and Scutaru to generate the Dekker master 
equation, we will operate under the assumption that this is also an appropriate form for 
a completely positive dynamical semigroup for unbounded operators. We will focus on the 
additional part of the master equation (Langevin equation) associated with dissipation, AL 
(AT). 

In Sec. y we will examine the classical Langevin equations for a particle coupled to a 
harmonic oscillator bath, which will help clarify the connection between first principle clas- 
sical and quantum calculations, as well as the issues which arise for quantum generalizations 



of classical systems. In Sec. |T|, we will construct AL by choosing a particular form for the 
operators {V^}, and establish some constraints on these operators from the requirement of 
linear dissipation. We will choose particular form for the operators which satisfy the con- 
straints, yield a master equation which is translationally invariant, and which correspond 
to an isotropic environment. The results are presented in both the Schrodinger and Heisen- 



berg forms. In Sec. [TV] we will derive and discuss fluctuation-dissipation relations for these 
models, as well as the time dependence of expectations of operators which are functions of 
position or momentum only. In Sec. [V], we compare our results for one dimension to the 



Dekker and Caldeira- Leggett master equations. We will comment on our results in Sec. VI 



II. CLASSICAL LANGEVIN EQUATIONS 

In this section, we review some aspects of Langevin equations which arise from interaction 
with a harmonic oscillator bath, essentially following Zwanzig ||. The composite system 
has a lagrangian given by 

L = -mx 2 — U (x) 

+ E^K-^-a,(x)) 2 ] (2.1) 

where x is the position of the system of interest, and are the coordinates of the oscillators 
comprising the bath. The interaction terms are of the form q ti m fJ- u 2 ,a fl (x), and the additional 
terms of the form ^m fJi u 2 a 2 l ('x) are associated with regularization of the potential. The 
equation of motion for an individual oscillator is given by 

% = -ufa + wJa M (x), (2.2) 

which is simply a driven simple harmonic oscillator. The formal solution to the equation of 
motion is 

%{t) = <? M (°) cosw M t + sin u^t 

+ / sinu; M (t - s)u; M a M (x(s))ds (2.3) 

J 



This becomes, after an integration by parts, 



qn(t) = qpWcosUpt+^—stiiUfJ + a^xfy) 



- coso; M (t - s)a fl (x(0)) 

<9a M (x(s))^ 



-J cos^(t- s) ^ x " Xjds. (2.4) 



The equation of motion for the particle is given (for the component X{) by 

Ht) = ~f| + E - a,(x))^. (2.5) 



Direct substitution of equation Eq. Q2.4D into Eq. ( J2.5[ ) yields an equation for motion for x 
which can be written 

d[/(x) ^ . , 
x * (*) = ^ + /i(x,t) 

— / r/jj(x(t), x(s); t — s)xj(s)ds, (2.6) 

« 

(summation over j implied) where f(x, t) is interpreted as a fluctuating force given by 

/i(x,t) =S{ m M a; J[(9A»( ) - a^(x(0))cosa; M t 

'^Bh^ y 11 }. (2.7) 



and the nonlinear dissipation kernel is given by 
?7y(x(t),x(s);t- s) = 

E"v? 1 ""^ - s) - (2 - 8) 

The initial probability distribution of the oscillators determines the statistics of the 
fluctuating force. The oscillators are assumed to be independent and in a thermal state, 
with 

(q,(0) - a„(x(0))) = (^(0)> = 0, (2.9) 

and 
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(-m^(g M (0) - « M (x(0))) 2 ) 



1 



(fe(o)) 2 ) 




(2.10) 



Under these conditions, we have the following fluctuation- dissipation relation: 



(f(x,t)> = 0, 



(/i( x >*)/i(y 5 s)) = ksTrjij^y; 



t-s). 



(2-11) 



With a typical Markov approximation, the width in time of r\ij is taken to be small, effectively 
a delta function. The fluctuating force is effectively white noise, and the effective Langevin 
equation can be written: 



The dependence of the dissipation on spatial correlations of the fluctuating forces has 
been lost. Since the particle can only be at one position at one time, these correlations 
play no role in the classical dynamics. Indeed, a wide variety of models can reduce to 
the same Langevin equation, if only in approximation. One could even assume a spatially 
homogeneous environment in which the rj^ is independent of x, and the resulting Langevin 
equation to get linear dissipation. A quantum mechanical description of the particle and 
its dynamics includes superpostions between different positions, resulting in interference 
between different particle trajectories. One effect of the interaction with the oscillator bath 
is that the noise will result in a loss of quantum coherence. This loss of coherence will 
depend upon the spatial correlations of the noise 0,0- l n the the following sections we 
will explore master equations which have linear dissipation and include spatial correlation 
effects, and thus are candidate models for the effects on a quantum system resulting from 
local coupling to a homogeneous and isotropic environment. 




(2.12) 



and the new fluctuation-dissipation relation is 




(2.13) 
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III. CONSTRUCTION OF THE MASTER EQUATION 



To construct the master equation, we focus on the non unitary part of the evolution, 
that is, 

The time dependence of expectation values can immediately be obtained via 

^ = ( ^ )+(i -[ ]>=(f) + (^,0]) + (Ar [0 l), (3.2) 

where (O) = Tr[pO]. We consider dissipation in the general case of d dimensions. For 
strictly linear dissipation, we require 

AL'\p] = -lj>, 
m 

AL*[x] = 0, (3.3) 

which will produce quantum Langevin equations corresponding to Eq. ( |1.1| ). One could also 
consider more general forms for the linear dissipation where rj is replaced with a dissipation 
tensor, and the classical frictional force components would be given by Fj = (summation 
over j), corresponding to an environment which was anisotropic. The anisotropic case could 
be addressed by a simple generalization of what we present here. 

The operators {V^} presumably could be expressed as functions of the position and 
momentum operators. This is the approach taken by Sandulescu and Scutaru citeScutaru, 
where {V^} were taken to be linear combinations of position and momentum. We will choose 
the particular form of {V^} to be at most linear in momentum, and write 

V, = A M (x) - B M (x) ■ p, (3.4) 

were {v4 M (x)} and {B M (x)} are as yet undetermined functions of x. The choice of this 
form is practically dictated by the requirement that the dissipation be linear, as any higher 
powers of momentum would appear as nonlinear orders of momentum in the first part of Eq. 
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(P-3|). We can now use Eq. (|3TT|) and Eq. (|3~3|) to establish some conditions on the functions 
{^(x)} and {5 M j(x)}. Evaluating Eq. (|3.1|) for a particular component of momentum, we 
take 

AL >«] = ^E^>n,^] + [V„, Pn ]V}, (3.5) 

substitute for {V^} with Eq. (|3.4|), and reorder the operators with all momentum terms to 
the right: 

[Pn] ~ 2, ( t {A * dx n dx n A »>> 2^ [ dx l {B - dx n dx n A » )] 
1 f dA^ dB^ dA\ dB^ t 

The resulting expression contains terms up to second order in momentum, and we equate 

each power to the corresponding term in Eq. ( pj.3[ ), which produces, after some simplification, 

the following conditions on {^(x)} and {B^x.)}: 

„ . dB uj 8B*. 
U '" dx n * 

i aA M db; oa; bb^ ^ 

2i ^ dx n dx n A » ^ dx n + dx n ^ ~ m m 
1 -dA u dA* d , , dA u dB* ui a ^ 

s " + * ^ (B < ~ <k a > )] = (37) 

To satisfy the second part of Eq. ( |3.3|) , we look at Eq. (|3.1|) for a particular component of 
the position operator to get 

^[x n ] = W[-^(Bt,B, n )] 

+7E[ fl J» fl Mj-^/«)Pi]}- (3-8) 

1 M 



The resulting constraint equations are 

Et^n^ - ^; y ^n)] = 0. (3.9) 

The conditions above by no means determine the functions {A M (x)} and {-^(x)}, but 
are constraints to be satisfied. The master equation we are generating is presumably a 
limiting form of the reduced dynamics of a composite quantum system, so we select a 
particular form of {A M (x)} and {5 M j(x)} to reflect the form of the interaction between the 
system of interest and the environment degrees of freedom. For a local interaction with a 
field, the environment degrees of freedom are the particular modes of the field. A reasonable 
choice for these functions are then 

{A,} = {a(k)e^} 

{B M } = {We**}, (3.10) 

where the index /i is replaced by a continuous d dimensional wave vector k, and sums over fj, 
are replaced by integrations over k. Some simplification occurs immediately since, with this 
choice, terms such as the product A* m B mi are independent of x. The conditions expressed 
in Eq. ( |3.7| ) and Eq. ( |3.9| ) reduce to 

J d d kb*{k)b j {k)k n = 

/ d d k2Re(a\k)b l (k))k n = 6 l A 
J m 

J d d k\a{k)\ 2 k n = 
J d d klm(a (k)fti(k)) = 

J d d jfclm(6*.(k)6 i (k)) = 0. (3.11) 

To satisfy these conditions, we need a(k) to be an even function of the components of k, 
and 6j(k) to be an odd function of hi, but an even function of the other components of k. 
To insure isotropy, we choose the form of the functions a(k) and o(k) as 
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a(k) = a(k) 

6,(k) = (3.12) 

These functions can then be seen to satisfy the conditions (Eqs. ( |3.7|) and (|3.7|)) by con- 
sidering the symmetry of the even limits for the integrals over particular components of k, 
along with the even or odd integrands. The dampening constant is given for this final choice 
by 

— = / d d k2Re(a(k)(3(k))kl 
m J 

r k 2 
= J d d k2Re(a(k)(3(k)) — . (3.13) 

We can now write the nonunitary part of the evolution equations for the particular choice 
of {V^}. For the Heisenberg picture (after some simplification using the symmetries of the 
functions a and (3) we have for an arbitrary operator O, 

AL* [O] = [ A^^(e-* k - x C>e ik - x - O) 



h 

+ J d d k^^-(k ■ pe^O^k ■ p - i{0, (k ■ p) 2 }) 
^ RB(q(fc)V(fc)) ({e ^ Qefa)k p}) 

^ flm("WV(fc)) ([e - fl c-x QeflMt>k , p]); (3 . 14) 

while for the Schrodinger picture we have 

AL[p] = - J d d k^^{p - e tk ->- ikx ) 

_ / d d k^^ (-{(k • p) 2 , p} - k ■ pe ikx pe- ikx k • p) 
J n 2 

_ J ^Re(a(fc)V(fc)) ({k p e , k . Xpe - Jk . x}) 

_ | ^ ^Im(a(fc)>(fc)) ([k p e , k . Xpe _, k . x]) _ (3 15) 

This is the main result of this paper: we have illustrated a method of constructing quantum 
master equations which are ohmic and translationally invariant. The translational invariance 
can easily be seen by replacing x with x' + r, where r is a c-number (p is unchanged by 
translation). The form of the nonunitary terms is unchanged. 
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IV. DISSIPATION, FLUCTUATION 
AND EXPECTATIONS 



In this section, we will examine the implications of the dissipative part of the evolution 
given by Eq. ( [3.14|) or Eq. (|3.15| ). by making use of the Heisenberg evolution equations, 
where the time dependence of expectation values is given by: 

Much immediate simplification can occur when one notes that e ik x (9 e -* k x j s simply a uni- 
tary transformation, corresponding to a translation of Kk in momentum space. Additional 
simplifications occur with the symmetries of the integrals over k and the particular form of 
a and as function of k = |k|. From this, one can readily verify that AL (p) = — -^-p. 

We now wish to look at the effect of the non unitary part of the evolution on kinetic 
energy. We find, with some simplification which arises from the isotropy of the model, that 

l 2m s J 2m 

H f^ k ml h k'-2^. (4.2) 
J a m 2m 

This can be interpreted as a relaxation to thermal equilibrium, under certain circumstances. 

We can make the identifications 

7S2 jl_ 



m J d 

-k B T = - [ d d k^^hk 2 , (4.3) 
2 7 J 2m v ' 

which determine an effective equilibrium temperature, provided 7 > 0. Since AL* is linear 

and since AL*(c — number) = 0, we can now rewrite Eq. (|4.2|) as 



A£ 'K£-5*» r >] = -*£-5*» r >- (4-4) 

The first term on the RHS of Eq. ( |3.15| ) can be identified as the effect of a random 
gaussian fluctuating potential. For a particle in a random potential with a translationally 
invariant two point correlation function with short relaxation time, approximated by 
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(V(x,t)V(y,s))=g(x-y)6(t-s), (4.5) 

the master equation gains an additional term which is given in the position representation 
by Hi 

(x|AL(p)|y) = -^G?(0) - <?(x - y))p(x, y). (4.6) 



The position representation of the first term on the RHS of Eq. (|3.15 ) is then given by 

\a(k)\ 2 



(x|/A^(p-e^>-' k -)|y) 



d d k H^l {1 _ e iH^y) )p ^ y); (4.7) 

so that after some simplification, we can identify this as the effect of a random gaussian 
potential with correlation function given by 

(V(x,t)V(y,s)) = 

h J d d k\a{k)\ 2 cos{k-{x-y))5(t-s). (4.8) 
In terms of fluctuating forces, we can write 

(f (x, t) • f (y, s)) = V, • V y <V(x, t)V(y, s)) 

= hj d d k\a{k)\ 2 k 2 cos{k-(x-y))5(t-s). (4.9) 

The fluctuation-dissipation relation can now be written using the identifications made in 
Eq. 



(f(x,t).f(x,s)) = (f(0,t)-f(0,s)> 

= 2m^k B T6(t- s). (4.10) 

We now wish to consider the expectations of more general quantities. For an arbitrary 
function of the position operator g(x), it is straightforward to examine the time dependence 
of its expectation value via Eq. Q4.1| ). Since x is unaffected by translations in momentum 
space, evaluating AL*[g(x)] via Eq. (|3.14|) becomes 
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AL*[g(x)] = 

+ / ^^-(k ■ PS (x)k ■ p - i{j(x), (k ■ p) 2 }) 
-/A Re '°'f^» ( {9 (x),k.p}) 

-/jfc <to(tt( f /?W) (b(x),k.ri), (4.11) 
which is further simplified via the symmetries of the integrals over k to 

AL>(x)] = / A^-^V^x) 

= § V 2 <?(x), (4.12) 

effectively defining the diffusion constant D. The spreading of the density operator can then 
be characterized by 

' ((X ^ Xo)2) = (^[^(x-x ) 2 ]) + D. (4.13) 
For an arbitrary function of momentum, h(p) we have 

e lkoc /i(p)e- lk - x = h(p + hk), (4.14) 

and AL*[h(p)] becomes 

AL>(p)] = 

/ d d k^^(h{p + hk)-h(p)) 

+ / d d k^^[{k • p) 2 (MP + hk) - h(p))\ 

- y <rfc v y v 77 (k • p)/i(p + ftk). (4.15) 



V. COMPARISON WITH DEKKER AND 
CALDEIRA-LEGGETT MASTER EQUATIONS 



It is useful to take a look at Sandulescu and Scutaru's construction of the Dekker master 
equation in order to make comparisons with our results. In section two, we were largely 
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followed their construction. Their choice for the operators {V^} were linear combinations of 
position and momentum operators, 

Vt = dix + kp, (5.1) 

while in our d-dimensional generalization, we considered operators of the form: 

Vk = a(k)e* koc - /3(k)e ik - x k • p. (5.2) 

It is possible to consider the Dekker equation as a low length scale limit of our result. One can 
imagine that this will correspond to a system which is well localized, perhaps by the potential 
of the system of interest. We can take the origin to be centered in the region of localization, 
noting that translations will not affect the form of AL. The length scale is determined 
the physical extent of this region of localization. A short length scale approximation can 
be considered as a long wavelength (for the environment modes) approximation via the 
"typical" wave numbers in Eq. ( |5.2|) . To make this approximation we expand in powers of 
k. To second order, we have 

e ik.x 0e -ik.x ^ o + i[k . x , O] - -[k • x, [k • x, O}}. (5.3) 
With simplifications (due to the symmetry in our choice of a and /3), AL becomes 

AL[p] = -J A^[k.x, [k-x,p]] 
-/^^[k-p,[k.p,p]] 
-Jd d k iR < a{ f m \ k.^{k. p , P }] 

In 1 dimension, this corresponds to the Dekker master equation, with 

AL[p] = -i^[x,{p,p}} + ^(D xp + D px )[x, [p,p\] 

DxX \p, \P,P\] ~ -np[z> [x,p]], (5.5) 
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where the parameters are given by 



2 
A 

2 

v(kYB(k))k 2 = 

2m 



X = J dkiRe{a{k)' (3{k))k 2 
(D xp + D px ) = J dkMm(a(k)*P(k))k 2 . (5.6) 



For the Caldeira-Leggett master equation, we have 

which can be viewed as a low momentum approximation (ignoring second order momentum 
terms in AL) of the Dekker master equation with D xp = D px = 0. If we make similar 
approximations for our master equation, without making short length scale approximations, 
we have 

AL[p] = - J d d k^^(p - e ikx pe- ikx ) 
_ J ^ Re(c(fc)>(fc)) ({k p ^ p ^ }) (5 g) 

For this case, the fluctuation-dissipation relations are of a more familiar form, relating the 
effective temperature to the correlation of the fluctuating potentials, where the equivalent 
relations to Eqs. ( |4.3| ) and( ^.10|) become 



and 



7 = 2— 

m 

d kB T^lf d ^^pl M ^ (5.9) 
2 7 J 2m 



(f(x,t).f(x,s)) = (f(0,*).f(0,s)) 

= 2m / y—kBTS(t — s) 

= Ar 1 h B T5{t-s). (5.10) 
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VI. DISCUSSION 



We have introduced a method of construction of models for linear (i.e. ohmic) dissipa- 
tion for a quantum system for arbitrary dimension. The model is completely specified by 
identifying the complex functions a(k) and b(k). Although we have focused on isotropic 
dissipation, anisotropic dissipation can readily be introduced by retaining the direction de- 
pendence of a and b on k. We have identified the phenomeno logical dissipation constant, and 
effective temperature, as well as the spatial correlation function for the fluctuating forces. 
Generalizations of this approach might also include considering other sets of non-trivial 
functions for {^(x)} and {13^ (x)}, so that the conditions expressed in Eqs. (|3.7|) and ( p.9|) 
are satisfied when averaged over some characteristic length scale. 

The models we have introduced are natural generalizations of the Dekker master equa- 
tion, and in one dimension, reduce to it in a short length scale approximation. Conceivable, 
one could construct a more general master equation using Dekker's noise operators by in- 
cluding spatial correlation effects. The Caldeira- Leggett master equation can be viewed as 
a short-length scale, small momentum approximation to our result. 

Caldeira and Leggett have pointed out that the only coupling to an oscillator bath 
which exactly produces ohmic dissipation is linear coupling (in both system and environment 
coordinates). However, one might expect that nonlinear coupling may provide dissipation 
which is linear under a wide range of approximations. The models we have proposed here 
are intended to anticipate the form of the dynamics which would result from local coupling 
to an environment. We have previously argued that the effect on the dynamics of quantum 
coherence may differ considerably from that which can be obtained by linear coupling. ||[7| 
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